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Abstract 



We present preliminary results of several analyses searching for the effects of 
CP violation and mixing in the decay of D° mesons. We find no evidence of CP 
asymmetry in five different two-body decay modes of the D° to pairs of light 
pseudo-scalar mesons: Acp{K+K~) = (+0.05±2.18±0.84)%, A C p(tt + tt-) = 
(+2.0±3.2±0.8)%, A C p(K^°) = (+0.1±1.3)%, A C p(tt°tt ) = (+0.1±4.8)% 
and Acp(K^Kg) = (—23 ± 19)%. We present the first measurement of the 
rate of wrong-sign D° K+ir~ir° decay: R ws = 0.0043ig;g{£J ± 0.0007. 
Finally, we describe a measurement of the mixing parameter yep = ^ 
searching for a lifetime difference between the CP neutral K + tt~ final state 
and the CP even K + K~ and 7r + 7r~ final states. Under the assumption that 
CP is conserved we find y C p = -0.011 ± 0.025 ± 0.014. 
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I. INTRODUCTION AND MOTIVATION 



The study of mixing in the K° and B® sectors has provided a wealth of information 
to guide the form and content of the Standard Model. In the framework of the Standard 
Model, mixing in the charm meson sector is predicted to be small |l| , making this an excellent 
place to search for non-Standard Model effects. Similarly, measurable CP violation (CPV) 
phenomena in strange 0|J and beauty [|]-|(|] mesons are the impetus for many current and 
future experiments [|7|-|1~0|1. The Standard Model predictions for CPV for charm mesons are of 



the order of 0.1% (TT|, with one recent conjecture of nearly 1% [12]. Observation of CPV in 



charm mesons exceeding the percent level would be strong evidence for non-Standard Model 
processes. 

A D° can evolve into a D° through "ordinary" on-shell intermediate states, or through 
off-shell intermediate states such as as those that might be present due to new physics. We 
denote the amplitude through the former (latter) states by — iy (x), in units of T D o/2 [[13 



The Standard Model contributions to x are suppressed to \x\ ~ tan 6 C ~ 5% and the 
Glashow-Illiopolous-Maiani fl4| cancelation could further suppress \x\ down to 10~ 6 — 10~ 2 . 
Many non-Standard Model processes could lead to \x\ > 1%. Contributions to x at this level 
could result from the presence of new particles with masses as high as 100 - 1000 TeV ||15|| . 
Signatures of new physics include \x\ ^> \y\ and CP violating interference between x and y 
or between x and a direct decay amplitude. 

Wrong sign (WS) processes, such as D° — > K + tt^tt°, can proceed directly through doubly 
Cabbibo-suppressed decay (DCSD) or through mixing and subsequent Cabbibo favored decay 
(CFD). Both DCSD and mixing followed by CFD can contribute to the time integrated 
WS rate, -Rws — (r + f)/2 and the inclusive CP asymmetry A = (r — f)/(r + f), where 
r = T (D° — > /) /r (^D° —►/), t is the charge conjugated quantity, and / is a wrong-sign 
final state, such as fT + 7r~7r°. 

The different contributions to i?ws and A can be separated by studying the proper decay 
time dependence of WS final states, as we have done in D° — > K + n~ fL6| , and has been done 
by FOCUS in |TTJ. The differential WS rate relative to the right-sign (RS) process, in time 
units of the mean D° lifetime, t D o = (415 ± 4)fs |L8|, is r{t) = [R D + y/lfoy't + l/4(x' 2 + 
l/' 2 )t 2 ]e _ * Ul9| , where y' = y cos 5 — x sin 5 and x' = x cos 5 + y sin 5, where Rn is the relative 



rate of DCSD and S is the strong phase between the DCSD and CFD amplitudes. There 
are theoretical arguments that 5 should be small ||20|| , although this should not be taken for 
granted. The coefficient of the term quadratic in t is proportional to the relative rate of 
mixing, R M = l/2(x 2 + y 2 ) = l/2(x' 2 + y' 2 ). 



II. GENERAL EXPERIMENTAL METHOD 

All of the analyses discussed herein, unless otherwise stated, use the same data set and 
reconstruction techniques described below. The data set was accumulated between February 
1996 and February 1999 and corresponds to 9.0 fb _1 of e + e~ collision data at ^fs ~ 10.6 
GeV provided by the Cornell Electron Storage Ring (CESR). The data were recorded by the 
CLEO II detector |2IJ upgraded with the installation of a silicon vertex detector (SVX) |^2[ 
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and by changing the drift chamber gas from an argon-ethane mixture to a helium-propane 
mixture |23|. The upgraded configuration is referred to as CLEO II. V. 

The Monte Carlo simulation of the CLEO II. V detector was based upon GEANT |[24|| , 
and simulated events were processed in the same manner as the data. 

The D° candidates are reconstructed through the decay sequence D* + — > D°7i+ |£5[] . The 
charge of the slow pion (71"+) tags the flavor of the D° candidate at production. The charged 
daughters of the D° are required to leave hits in the SVX and these tracks are constrained 
to come from a common vertex in three dimensions. The trajectory of the D° is projected 
back to its intersection with the CESR luminous region to obtain the D° production point. 
The 7r+ is refit with the requirement that it come from the D° production point, and the 
confidence level of the x 2 °f this refit is used to reject background. 

The energy release in the D* — > D°n^ decay, Q = M* — M — m n , obtained from the above 



technique is observed to have a width of oq = 190 ± 2 keV ||26|| , which is a combination of 
the intrinsic width and our resolution, where M and M* are the reconstructed masses of the 
D° and D* + candidates respectively, and is the charged pion mass. The reconstruction 
technique discussed above has also been used by CLEO to measure the D* + intrinsic width, 
r D *+ = 96 ± 4 ± 22 keV (preliminary) ||27|| . In the mixing analyses described below, the 
distribution of candidates in the Q vs M plane are fit to determine both right-sign and 
wrong-sign yields. 

We calculate t using only the vertical component of the D° candidate flight distance. This 
is effective because the vertical extent of the CESR luminous region has cr vcrtical = 7fim P| . 
The resolution on the D° decay point (x v , y v , z v ) is typically 40/mi in each dimension. We 
measure the centroid of the CESR luminous region (xb, yb, Zb) using e + e~ — > qq (q = udscb) 
events from sets of data with integrated luminosity of several pb _1 , obtaining a resolution 
on the centroid of less than 5/zm. We express t as t = M/p vert icai x (y v — yb)/ (cr D a), where 
^Vertical is the vertical component of the total momentum of the D° candidate. The error in 
t, a t , is typically 0.4 (in D° lifetimes), although when the D° direction is near the horizontal 
plane a t can be large. 



III. CP VIOLATION IN D° DECAY 

CP Violation in charm meson decay is expected to be small in the Standard Model, which 
makes charm meson decay a good place to look for non-Standard Model effects. Cabbibo 
suppressed charm meson decays have all the necessary ingredients for CP violation - multiple 
paths to the same final state and a weak phase. However, in order to get sizable CP violation, 
the final state interactions need to contribute non-trivial phase shifts between the amplitudes. 
Large final state interactions are a likely reason why the prediction for the ratio of branching 
ratios of (D° — > K + K~)/ (D° — > 7r + 7r~) yields a value roughly half of the observed value 
hence these may provide a good hunting ground for CP violation. 



Previous searches for mixing-induced [16] or direct [ 2£ , 18] CP violation in the neutral 



charm meson system have set limits of ~ 30% or a few percent, respectively. We present 
results of searches for direct CP violation in neutral charm meson decay to pairs of light 
pseudo-scalar mesons: K + K~, 7r + 7r~, /^s7t°, 7r°7r° and K^K^. 
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A. Search for CP violation in D° — ► K + K and D° — > 7r + 7r decay 
The asymmetry we want to measure, 

" r(£>°-/) + r(se-/) 

can be obtained from the asymmetry 

_ r p*+ -> 7T+/) - r (d*- -» 7T-/) 

r( jD ^^7r+/) + r( J D-^7r-/) 

The slow pion and D° are produced by the CP-conserving strong decay of the D* + , so the 
slow pion serves as an unbiased flavor tag of the D°. The decay asymmetry can be obtained 
from the apparent production asymmetry shown above because the production of D is 
CP-conserving. 

The asymmetry result is obtained by fitting the energy release (Q) spectrum of the 
D* + — > D°7T^ events. The D° mass spectra are fit as a check. The background-subtracted 
Q spectrum is fit with a signal shape obtained from K + n~ data and a background shape 
determined using Monte Carlo. Figure [1] shows our fits to both D — > 7r + 7r~ and D — > K + K~ 
distributions found in data, for both signs of the slow pion flavor tag. The parameters of the 
slow pion dominate the Q distribution, so all modes have the same shape. We do the fits 
in bins of D° momentum to eliminate any biases due to differences in the D° momentum 
spectra between the data and the MC. The preliminary results are A(K + K~) = 0.0005 ± 
0.0218(stat) ± 0.0084(syst) and A(it + it-) = 0.0195 ± 0.0322(stat) ± 0.0084(syst). 

We use many different variations of the fit shapes, both empirical and analytical, to assess 
the systematic uncertainties due to the fitting procedure (0.69%). We also consider biases due 
to the detector material (0.07%), the reconstruction software (0.48%), and forward-backward 
acceptance variations (cc pairs are not produced symmetrically in the forward/backward 
directions in e + e~ collisions at y/s ~ 10.6 GeV, and the collision point was not centered 
exactly in the middle of the detector) (0.014%). 

The measured asymmetries are consistent with zero, and no CP violation is seen. These 
results compared to previous measurements [ff§| , |50]| are shown in Figures ^ and 



B. Search for CP Violation in D° i^vr u , D° -» 7r°7r° and D° -> K£K% decay 

This analysis [^l|] differs from the other analyses presented in this paper in some of its 
reconstruction techniques and in the data set used. The 7r°7r° and K^ir final states do 
not provide sufficiently precise directional information about their parent D° to use the 
intersection of the D° projection and the CESR luminous region to refit the slow pion as 
described in the general experimental technique section. The K^Kg final state is treated the 
same for consistency. This analysis uses the data from both the CLEO II and CLEO II. V 
configurations of the detector, totaling 13.7 fb _1 of e + e~ collision data at y/s ~ 10.6 GeV. 

The .Kg and n° candidates are constructed using only good quality tracks and show- 
ers. The tracks (showers) whose combined invariant mass is close to the Kg (it ) mass 
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FIG. 1. Q distributions for D° — ► K + K~ (left) and D° —>■ tt + 7t~ (right) candidates separated 
into samples tagged as D* + (top) and D*~ (bottom). Shown immediately to the right of each Q 
distribution are the point by point residuals of the fit. In each case the data points are shown with 
error bars and the solid line represents our fitting results. 
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FIG. 2. A(K + K ) from this analysis compared to previous results. 
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Summary of A cp (7i7i) 
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FIG. 3. A(tt + it ) from this analysis compared to previous results. 



are kinematically constrained to the K% (tt°) mass, improving the D° mass resolution. 
The tracks used to form Kg candidates are required to satisfy criteria designed to reduce 
background from D° — > 7t + tt~X decays and combinatorics. Candidate events with recon- 
structed D° masses close to the known D° mass are selected to determine the asymmetry, 
A(f) = [t(D° -»■ /) - T(W -> /)] / [V(D° -»■ /) + T(D° -> /)] . The Q distributions of the 
candidates in the three decay modes are shown in Figures f|, y, and A prominent peak 
indicative of D* + — > D°ir^ decay is observed in all three distributions. The total num- 
ber of D° and D° candidates for a given final state is determined as follows. We fit the 
Q distribution outside of the signal region and interpolate the fit under the signal peak to 
determine the background in the signal region. We subtract the background in the signal 
region from the total number of events there to determine the total number of signal events. 
The background shape is approximated as a non-relativistic threshold function with first and 
second order relativistic corrections B(Q) = aQ l l 2 + bQ 3 ^ 2 + cQ 5 ^ 2 . 

After background subtraction, we obtain 9099 ± 153 Kg7r° candidates, 810 ± 89 7r°7r° 
candidates, and 65 ± 14 K^K^ candidates. 

The difference in the number of D° and D° to a given final state is determined by taking 
the difference of the number of events in the signal region, and the asymmetry is obtained 
by dividing by the number of candidates determined above. This method of determining the 
asymmetry implicitly assumes that the background is symmetric. 

We have searched for any sources of false asymmetries: from the 7r+ finding (0.19%), 
from the fitting (0.5%), and from the backgrounds (0.35% i^vr , 0% tt°tt and 12% K%K$). 
We find no significant biases, but apply the measured corrections and add their uncertainties 
to the total. We obtain the results A(i^7T°) = (+0.1 ± 1.3)%, A(7r°7r°) = (+0.1 ± 4.8)% 
and A(K^Kg) = (—23 ± 19)% where the uncertainties contain the combined statistical 
and systematic uncertainties. All systematic uncertainties, except for the 0.5% uncertainty 
assigned for possible bias in the fitting method, are determined from data and would be 
reduced in future higher luminosity samples. 

All measured asymmetries are consistent with zero and no indication of significant CP vi- 
olation is observed. This measurement of A(Kgir ) is a significant improvement over previous 
results, and the other two asymmetries reported are first measurements. 



8 



3000 
2500 
2000 
1500 



S 500 
m 

in 
£1600 
a 

2 1400 
n1200 
1000 
800 
600 
400 
200 



(L 



_l_ 



(a) 



(b) 




0.005 0.010 0.015 0.020 0.025 

Q (GeV) 

FIG. 4. (a) Fitted Q distribution for D° — > K®ir candidates. The points with error bars are 
the data, the solid line represents the background and the dashed line shows the interpolation into 
the Q signal region, (b) The Q distributions for D° — ► K^ir (points) and D° — ► K^tt (histogram) 
candidates. 
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FIG. 5. (a) Fitted Q distribution for D° — > 7r°7r° candidates. The points with error bars are 
the data, the solid line represents the background and the dashed line shows the interpolation into 
the Q signal region, (b) The Q distributions for D — > (points) and D° — ► (histogram) 
candidates. 
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IV. FIRST OBSERVATION OF WRONG-SIGN D 



o 



K~ 



TT TT 



DECAY 



D° — ► Ktttt candidates are reconstructed using the selection criteria described in Sec- 
tion [TT], with additional requirements specific to this analysis. In particular, tt° candidates 
with momenta greater than 340 MeV/c are reconstructed from pairs of photons detected in 
the Csl crystal calorimeter. Backgrounds are reduced by requiring specific ionization of the 
pion and kaon candidates to be consistent with their respective hypotheses. 

The RS mode was recently studied by CLEO |32| and found to have a rich Dalitz structure 
consisting ofp(770)+, K*(&92)~, K*(892)°, p(1700)+ ^o(1430)°, #o(1430)- and K*( 1680)- 
resonances and non-resonant contributions. Recent theoretical predictions based on U-spin 
symmetry arguments [R3 suggest that the wrong sign (WS) channel will have a different 
resonant substructure than the right sign (RS). This analysis allows for different average 
WS and RS efficiencies resulting from different WS and RS kinematic distributions in the 
calculation of the WS rate: 



ttws — 



V 



WS 



RS 



The ratio of yields in Eq. (j^) is measured by performing a maximum likelihood fit to the 
two-dimensional distribution in m(KTciT) and Q. The signal distribution in these variables 
is taken from the RS data. The backgrounds are broken down into three categories: 1) RS 
D° — > Knir decay combined with an uncorrelated its, 2) combinations from e + e~ uu, 
dd, and ss, and 3) combinations from charm particle decays other than the RS or WS signal 
modes. The background distributions are determined using a large Monte Carlo sample, 
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which corresponds to approximately eight times the integrated luminosity of the data sample. 
The Q-m(K7nr°) fit yields a WS signal of 38 ± 9 events and a ratio N WS /N RS = 0.0043±^. 
Projections of the data and fit results in slices through the signal region in each variable 
are shown in Fig. [7[ The statistical significance of this signal is found to be 4.9 standard 
deviations. 

The average efficiency ratio in Eq. ([!]) is determined using a fit to the Dalitz plot variables 
m 2 (K + TT~) and m 2 (K + TT°) in the WS data. In this fit, the amplitudes and phases are initial- 
ized to the RS values, and those corresponding to the K*(892) + and K*(892)° resonances 
are floated relative to the dominant p(770)~ and other minor contributions. Combining the 
square of the fitted amplitude function with a parameterization of the efficiency determined 
using a large non-resonant Monte Carlo sample, we measure an average efficiency ratio of 
1.00 ± 0.02(stat). Studies are under way to examine the extent and the significance of this 
surprising similarity between the RS and WS Dalitz plots. 

The dominant systematic errors in this analysis come from the uncertainty in the Monte 
Carlo background distributions (14%), uncertainty in the amplitudes and phases that are 
fixed in the fit (8%), and uncertainty in the background Dalitz plot distributions (3%). 

Several powerful checks of the Q-m(iC7T7r ) and Dalitz plot fits are performed in order to 
verify the validity of these results. Fits using specific background regions of the Q-m(PC7T7r ) 
plane test the sensitivity to the Monte Carlo background distributions. The WS Dalitz 
plot is also fit using hypotheses which include one of the iT > (892) + , i^*(892)°, or p(770)~ 
resonances, in order to provide an upper limit on the error due to this fit. Only the i^*(892)° 
hypothesis leads to an efficiency ratio that differs significantly from one, but this hypothesis 
is strongly disfavored by the data. 

We measure the wrong sign rate to be 

R ws = 0.0043^;^ (stat) ± 0.0007 (syst) (preliminary). (2) 

This result is consistent with the CLEO II. V [[llj and FOCUS jT7| D° — > Kn measurements. 
This measurement of Rws can be used to obtain limits on Rdcsd as a function of y' within 
the limits on x' set by CLEO, as shown in Figure [8]. Work is in progress to use the lifetime 
distribution of this sample to yield independent limits on x', y', and Rdcsd- 



V. SEARCH FOR CP DEPENDENT LIFETIME DIFFERENCES DUE TO D° - D° 

MIXING 



In the limit of no CP violation in the neutral D system we can write the time dependent 
rate for D — > /, where / is a CP eigenstate, as 

R(t) oc e - tr{1 ~ vcpricp) (3) 

where T is the average D width, t]cp is the CP eigenvalue for /, and 

Vcp = 2p ( 4 ) 

where Ar is the width difference between the physical eigenstates of the neutral D. In the 
limit that CP is conserved in charm decays, yep = y, where y is the mixing parameter defined 
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FIG. 7. Results of the fit to the data to determine Nws/Nrs- Projections in the variables a) 
Q and b) m(KTnr ), after selecting the signal region (within two standard deviations) in the other 
variable. 
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candidates. The curve is the fit described in the text. 
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TT 



in Section |. With the CP asymmetry for KK and nn consistent with zero as discussed in 



Section |TJ, this limit is well motivated and we will simply use y for the rest of this section. 
We can then express y as 



y 



'CP 
TCP+ 



(5) 



where v^jp is the lifetime of a CP neutral state, such as Kn, and tcp+ is the lifetime of a 
CP even state, such as KK and tttt. Thus to measure y we simply take the ratio of the 
lifetimes of D° — ► Kir to D° — > KK and -kti. Since the final states are very similar, our 
backgrounds are small and cross-feed among the final states is negligible many of the sources 
of uncertainty cancel in the ratio. A similar analysis has recently been published by FOCUS 
comparing the Kn and KK final states ||35|| . 

The left-most plots in Figures |^, (TTj, and [TT] show the candidate D° masses for the three 
samples selected as described in Section [TT]. Selecting events in a narrow region around the 
D mass, we fit their proper time distributions with an unbinned maximum likelihood fit. 
The proper time is measured as described in Section [TT]. The signal likelihood function is an 
exponential convolved with three resolution Gaussians. The width of the primary Gaussian 
is due to the propagation of errors from the track fit to the flight distance and momentum 
for each candidate. The second and third Gaussians represent candidates that have been 
mismeasured by the addition of spurious tracking hits or due to hard non-Gaussian scatters 
in the material of the detector. The first of these has its width determined in the fit to the 
copious Kn sample while the second is fixed to a large value, 8 ps. The relative contribution 
of the two mismeasured signal resolutions is determined using the Kn data sample and 
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ii—i. 



1.825 1.845 1.865 1.885 1.905 -4 4 8 

7171 Mass (GeV/c") 7C7I Proper Time (ps) 

FIG. 11. The invariant mass (left) and proper time (right) distributions for D° — » ?r 
candidates. The curve is the fit described in the text. 
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TABLE I. Summary of the lifetime fits. The parameters are those described in the text. Note 
that we have constrained the candidates to a D° mass of 1.86514 GeV, the Monte Carlo corrected 
weighted average of the data. The signal lifetime is highly dependent on the D° mass used in the 
constraint, while the lifetime difference is not. This technique yields the smallest uncertainty in y, 
but it not optimal for measuring the absolute D° lifetime and was not used in [34]. 



Parameter 


Kit 


KK 


7T7T 


Number of Signal 


20272 ± 178 


2463 ± 65 


930 ± 37 


Tsig (PS) 


0.4046 ± 0.0036 


0.411 ±0.012 


0.401 ±0.017 


Background Frac (%) 


8.8 ±0.2 


50.7 ±0.7 


29.1 ± 1.3 


Background Life Frac (%) 


81.0 ±4.8 


85.7 ±2.9 


32.2 ±7.5 


Uack (PS) 


0.376 ± 0.030 


0.436 ± 0.020 


0.56 ±0.15 


fmis % 


3.8 ±0.9 


Fixed 


Fixed 


<?mis (PS) 


0.590 ± 0.079 


Fixed 


Fixed 



fixed for the KK and rcn samples. According to our simulation, the fraction of the second 
Gaussian is about 4% of the well measured signal and the third Gaussian is less than 0.1% 
of the signal. The probability for a candidate to be signal is determined by its measured 
mass, and is based on the fit to the mass distributions. 

Background is considered to have contributions with both zero and non-zero lifetimes. 
All parameters that describe the background are allowed to float in the fits except for the 
width of the widest Gaussian which is fixed to 8 ps. 

The right-most plots in Figures [| [H], and [II] show the lifetime fits to the three samples 
and the results are summarized in Table |. 

Note that the background in all the samples has a large component that has a lifetime 
consistent with the D° lifetime. This agrees with the prediction of our simulation that the 
background with lifetime is dominated by misreconstructed fragments of charm decays. 

We calculate y separately for the KK and tttc samples. Systematic uncertainties are 
dominated by the statistical uncertainty in a Monte Carlo study used to determine small 
corrections, consistent with zero, that are applied to the measured result to account for 
differences between measured and generated values of the lifetimes (±0.009). Additional 
significant systematic uncertainties come from variations in description of the background 
(±0.008), uncertainties in our model of the proper time resolution (±0.005), and details of 
the fit procedure (±0.005), where the listed values are the contribution to the average result. 
Our preliminary results are 

y KK = -0.019 ± 0.029(stat) ± 0.016(syst) (6) 

and 

y n7T = 0.005 ± 0.043(stat) ± 0.018(syst). (7) 
We form a weighted average of the two to get 

y = -0.011 ± 0.025(stat) ± 0.014(syst) (preliminary). (8) 
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D°-D° Mixing Limits 



o.i 



i CLEO D°->KV 
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0.1 0.2 

x' (r D /2) 

FIG. 12. Allowed regions, at 95% CL, in the x' vs y' plane for some recent results. 



This result is consistent with zero. It is also consistent with our previous result from 
wrong sign Kit |16[], the FOCUS results from wrong sign Kit and the FOCUS results 



from Kir versus KK |35] . This result is shown along with some of the previous measurements 
in Figure [T^. 



VI. SUMMARY 

We present preliminary results of several analyses searching for the effects of CP violation 
and mixing in the decay of D° mesons. We find no evidence of CP asymmetry in five differ- 
ent two-body decay modes of the D° to pairs of light pseudo-scalar mesons: A C p(K + K~) = 
(0.05 ± 2.18 ± 0.84)%, A CP (n+n-) = (2.0 ± 3.2 ± 0.8)%, A cp (K°ti ) = (+0.1 ± 1.3)%, 
A C p(n 7r°) = (+0.1 ± 4.8)% and A CP (KgK%) = (-23 ± 19)%. We present the first mea- 
surement of the rate of wrong-sign D° -> K + n-n° decay: R ws = 0.0043lg:8oio ± 0.0007. 
Finally, we describe a measurement of the mixing parameter y = ^ by searching for a 
lifetime difference between the CP neutral K + ir~ and the CP even K + K~ and 7r + 7r~ final 
states. Under the assumption that CP is conserved we find y = —0.011 ± 0.025 ± 0.014. 
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